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Abstract 

The tensor electric polarizability of the deuteron gives important information about spin- 
dependent nuclear forces. If a resonant horizontal electric field acts on a deuteron beam circulating 
into a storage ring, the tensor electric polarizability stimulates the buildup of the vertical polar- 
ization of the deuteron (the Baryshevsky effect). General formulas describing this effect have been 
derived. Calculated formulas agree with the earlier obtained results. The problem of the influ- 
ence of tensor electric polarizability on spin dynamics in such a deuteron electric-dipole-moment 
experiment in storage rings has been investigated. Doubling the resonant frequency used in this 
experiment dramatically amplifies the Baryshevsky effect and provides the opportunity to make 
high-precision measurements of the deuteron's tensor electric polarizability. 
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I. INTRODUCTION 



Electric and magnetic polarizabilities are important properties of deuteron and other 
nuclei. Tensor electric and magnetic polarizabilities are defined by spin interactions of 
nucleons. In particular, measurement of the tensor electric polarizability of deuteron gives 
an important information about an interaction between spins of nucleons and provides a 
good possibility to examine the theory of spin-dependent nuclear forces. 

The method of determination of this important electromagnetic property of deuteron has 
been proposed by V. Baryshevsky et al. [1, 2, 3]. If an electric field acts on a deuteron beam 
circulating into a storage ring, the presence of the tensor electric polarizability leads to the 
appearance of an interaction quadratic in spin. When the electric field in the particle's 
rest frame oscillates at the resonant frequency, the effect similar to the nuclear magnetic 
resonance (NMR) takes place. This effect stimulates the buildup of the vertical polarization 
(BVP) of deuteron beam [1, 2, 3]. 

In the present work, we derive general formulae describing the BVP caused by the tensor 
electric polarizability of deuteron in storage rings (the Baryshevsky effect). Another effect 
defined by the tensor magnetic polarizability of deuteron consists in the spin rotation in 
the horizontal plane at two frequencies instead of expected rotation at the g— 2 frequency 
[1, 2, 3]. In the above cited works, the approach based on equations defining dynamics of 
polarization vector and polarization tensor has been used. To check obtained results and 
develop a more general theory, we follow the quite different method of spin amplitudes (see 
Refs. [4, 5]). In the present work, this method is partially changed. We use the matrix 
Hamiltonian for determining an evolution of spin wave function. 

The Baryshevsky effect should be taken into account at a search for the electric dipole 
moment (EDM) of deuteron [1, 2, 3]. An existence of the electric dipole moment of deuteron 
also leads to the BVP. It is planned to measure the BVP in the deuteron EDM experiment 
in storage rings [6, 7, 8]. Since the Baryshevsky effect can imitate the existence of the EDM, 
the spin dynamics caused by this effect should be investigated in detail. 

We use the relativistic system of units h = c = 1. 

In the next section, we review the main aspects of the used Hamiltonian approach in the 
method of spin amplitudes. In Section III, we briefly discuss the form of the Hamilton opera- 
tor in a cylindrical coordinate system. Section IV is devoted to the calculation of corrections 



2 



to the Hamilton operator for the tensor polarizabilities of deuteron. A solution of the matrix 
Hamilton equation is given in Section V. In Section VI, we calculate the spin dynamics ex- 
pressed by the evolution of the vertical component of polarization vector. A detailed analysis 
of new experimental possibilities for measuring the tensor electric polarizability of deuteron 
is performed in Section VII. Section VIII is dedicated to the differentiation of effects of the 
EDM and the tensor electric polarizability in the deuteron EDM experiment. Finally, in 
Section IX we discuss previously obtained formulae and summarize the main results of the 
work. 



II. HAMILTON! AN APPROACH IN THE METHOD OF SPIN AMPLITUDES 

The method of spin amplitudes uses quantum mechanics formalism to more easily describe 
spin dynamics (see Refs. [4, 5]). For spin-1/2 particles, it is mathematically advantageous to 
use this formalism because transporting the two-component spin wave function (spinor) \P 
is simpler than transporting the three-dimensional polarization vector P. The relationship 
between them is given by the expectation value of the Pauli spin vector <r: 

( 

P = ^cr^, * = 



(1) 



C+i/ 2 (t) 

\ c_i /2 (0 

where C+i/ 2 (i) and C_i/ 2 (i) are the time-dependent amplitudes. Together with the identity 
matrix the Pauli matrices generate an irreducible representation of the SU(2) group. 

Algebraically, the SU(2) group is a double cover of the three-dimensional rotation group 
SO (3). Therefore, the formalism based on the Pauli matrices is applicable to particles/nuclei 
with arbitrary spin if an effect of spin rotation is analyzed. The spin rotation can also be 
exhaustively described with the polarization vector P which is defined by 

<Sj> . . 

P%= — g — , hJ = x,y,z, (2) 

where Si are corresponding spin matrices and S is the spin quantum number. The polariza- 
tion vector being an average spin is a strictly classical quantity (see Ref. [9]) whose evolution 
can be investigated in the framework of classical spin physics. 

Particles with spin S > 1 also possess a tensor polarization. Main characteristics of such 
a polarization are specified by the polarization tensor which is given by [10] 
p _ 3<S i S j + S j S i >-2S(S+l)6 ij . . _ 

13 ~ 2S(2S-1) ' 1 >3- X >y> z - W 
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The polarization tensor satisfies the conditions Pij = Pji and P xx + P yy + P zz = 1 and 
therefore has five independent components. Additional tensors composed of products of 
three or more spin matrices are needed only for the exhaustive description of polarization 
of particles/nuclei with spin S > 3/2. 

The spin matrices for spin-1 particles have the form 

^ 1 \ ( -1 \ 



1 

71 



1 1 

V° 1 V 



So, 



V2 



1 

V° 1 



-i 
J 



s z 



/ 1 ^ 


V° 



(4) 



Possibly, the nontrivial spin dynamics predicted in Refs. [1, 2, 3] and conditioned by the 
tensor electric polarizability of deuteron is the first example of importance of spin tensor 
interactions in the physics of polarized beams. Tensor interactions of deuteron can also be 
described with the method of spin amplitudes. In this case, three-component spinors and 
3x3 matrices should be used. The method of spin amplitudes is mathematically advanta- 
geous because transporting the three-component spinor is much simpler than transporting 
the three-dimensional polarization vector P and five independent components of the polar- 
ization tensor P^ together. 

We follow the traditional quantum mechanical approach perfectly expounded by R. Feyn- 
man [11] and use the matrix Hamilton equation and the matrix Hamiltonian H for deter- 
mining an evolution of the spin wave function: 

' Cx(f) \ 
* = C (f) 

where H is 3 x 3 matrix, \1/ is the three-component spin wave function (spinor), 
and i,j = 1, 0, —1. In this equation, are matrix elements of the Hamilton operator Ti. 
A determination of spin dynamics can be divided into several stages, namely 

i) a solution of Hamilton equation (5) and a determination of eigenvalues and eigenvectors 
of the Hamilton matrix H; 

ii) a derivation of spin wave function consisting in a solution of a set of three linear 
algebraic equations; 

iii) a calculation of time evolution of polarization vector and polarization tensor. 



•IT = H% 



Hij =< i\H\j >, 



(5) 



4 



III. HAMILTON OPERATOR IN A CYLINDRICAL COORDINATE SYSTEM 



The spin dynamics can be analytically calculated when a storage ring is either circular or 
divided into circular sectors by empty spaces. In this case, the use of cylindrical coordinates 
can be very successful. The particle spin motion in storage rings is usually specified with 
respect to the particle trajectory. Main fields are commonly defined relatively the cylindrical 
coordinate axes. Therefore, the use of the cylindrical coordinates considerably simplifies an 
analysis of spin rotation in the horizontal plane (g— 2 precession) and other spin effects. 
Equation of spin motion in storage rings in a cylindrical coordinate system has the form [12] 



7=^x5, «. = -- \aB - -^T/3(/3 • B) + - a ) (f3 k E) 

at m I 7 + 1 \Y — 1 



1 

+ - 

7 



*||-^(0x£O| 



+Ue-^—(3((3-E)+(3xB 
2 V 7 + 1 



(6) 



where a — (g — 2)/2, g — 2fim/(eS), rj = 2dm/ (eS), and d is the EDM. The sign || means 
a horizontal projection for any vector. In this work, we do not consider effects caused by 
perturbations of particle trajectory investigated in Ref. [12]. The quantity u> a is also equal 
to [12] 

u a = n - 0e z , (7) 

where ft is the Thomas-Bargmann-Michel-Telegdi (T-BMT) frequency [13] corrected for the 
EDM [12, 14, 15, 16] and <fie z is the instantaneous angular frequency of orbital revolution. 

If we used the Hamiltonian of the particle with the EDM given in the lab frame [16], we 
would present the matrices S p and S<f, in the form 

S p = S x cos (j) + S y sin 0, — — S x sin + S y cos0. (8) 

However, this representation of spin matrices S p , S<p leads to cumbersome calculations be- 
cause the azimuth <f> defined by a particle position is time- dependent. Therefore, it is helpful 
to consider spin effects in the frame rotating at the instantaneous angular frequency of or- 
bital revolution which is almost equal to the cyclotron frequency. In this frame, the motion 
of particles is relatively slow because it can be caused only by beam oscillations and other 
deflections of particles from the ideal trajectory. The equation of spin motion in the rotating 
frame coincides with that in the cylindrical coordinate system because the horizontal axis 
of this system rotates at the instantaneous angular frequency of orbital revolution. 
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The Hamiltonians of the particle in the rotating frame and in the lab one (TL and TCi a b, 
respectively) are related by [18] 

n = n lab -s-uj, (9) 

where u) is the observer's proper frequency of rotation [19]. In the considered case, this 
frequency coincides with the instantaneous angular frequency of orbital revolution. The 
relation between the Hamiltonian in the lab frame and the T-BMT frequency corrected for 
the EDM is given by 

'Hiab = 7~Lo + S ■ O, 

where Ho is a sum of spin-independent operators. Therefore, the Hamiltonian in the rotating 
frame has the form 

H = H + S-uj a , (10) 

where ui a is defined by Eq. (6). Evidently, Hamiltonian (10) is consistent with Eq. (6). 

The particle in the rotating frame is localized and ideally is in rest. Therefore, we can 
direct the x- and y-ax.es in this frame along the radial and longitudinal axes, respectively. 
This procedure is commonly used (see Refs. [4, 5, 10]) and results in the direct substitution 
of spin matrices (4) for S p and Sf. 



S p — S x 



1 

71 



1 

1 1 
1 



\ 



Srf) — Sn 



I 

V2 



( 



0-10 
1 -1 
1 



\ 



(11) 



The matrix S z remains unchanged. Similar substitution can be performed for the matrices 
U ir 

The use of definition (11) strongly simplifies calculations. 

If an interaction causes a correction to Hamilton operator (10) and does not appreciably 
influence the particle motion, this correction is the same in the lab frame and the rotating 
one. 

In this work, we suppose B to be upward. For the deuteron, a < and {u) a ) z > 0. 
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IV. CORRECTIONS TO THE HAMILTON OPERATOR FOR TENSOR POLAR- 
IZABILITIES OF DEUTERON 

Corrections to the Hamilton operator for deuteron polarizabilities contain scalar and 
tensor parts. The scalar part is spin-independent and can be disregarded. General formulae 
used in Refs. [1, 2, 3] are within first-order terms in the normalized velocity (3. In the 
present work, we derive exact formulae for the configuration of main fields related to the 
resonant deuteron EDM experiment (see Refs. [6, 7]). Because the Lorentz factor is planned 
to be 7 = 1.28 [8] in this experiment, exact relativistic formulae are needed. 

Within first-order terms in (3, the interaction Hamiltonian depending on the electric and 
magnetic polarizabilities is given by 

V = V e + V m = -\a lk E[E' k - l -(3 ik B[B' k , 
E' = E + (3 x B, B' = B — (3 x E, 

where a^. and are the tensors of electric and magnetic polarizabilities, E' and B' are 
effective fields acting on a particle (fields in the particle's rest frame, i.e., in the rotating 
frame). In this approximation, the spin-dependent part of the Hamiltonian defined by the 
tensor electric and magnetic polarizabilities is equal to [1, 2, 3] 

V = -aT(S-E') 2 -p T (S-B') 2 , (13) 

where and (3? are the tensor electric and magnetic polarizabilities, respectively. 

The Baryshevsky effect takes place when one stimulates coherent longitudinal oscillations 
of particles at a resonance frequency. The angular frequency of forced longitudinal oscilla- 
tions, uj, is equal to the difference between two radio frequencies (see Ref. [8]). It should be 
very close to the angular frequency of spin rotation (g— 2 frequency), u , and close to the 
eigenfrequency of free synchrotron oscillations (synchrotron frequency) [8]. The resonant 
electric field in the particle's rest frame possesses the oscillating longitudinal component E^ 
defined by the Lorentz transformation of longitudinal electric field met by the oscillating 
particles and the radial one E' caused by the Lorentz transformation of vertical magnetic 
field. The latter component has a resonance part because of the modulation of the particle 
velocity. In the present work, we consider only effects of resonant fields on the BVP in 
ideal conditions and disregard systematical errors listed in Section VII. Thus, we take into 
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consideration the constant vertical magnetic field and the oscillating longitudinal electric 
one. 

The relativistic formulae for the fields in the particle's rest frame, E' and B', have the 
form 

E' = f3 1 B z e p + E^ B' = 1 B z e z . (14) 

The fields and the electromagnetic moments in the lab frame are unprimed. 

Induced electric and magnetic dipole moments in the particle's rest frame caused by the 
tensor polarizabilities are equal to 

d! = a T {S, (S • E')}, m> = (3 T {S, (S • B')}, (15) 

where {...,...} means an anticommutator. 

The correction to the Hamilton operator for the tensor polarizabilities of deuteron is 
equal to 

V = V lab = ~(d-E + m-B). (16) 

This correction does not change the angular frequency of orbital revolution. According to 
Eq. (9), the correction is the same in the rotating frame and the lab one. Since the induced 
dipole moments are proportional to the effective fields in the particle's rest frame, the factor 
1/2 appears. 

To obtain the dipole moments in the lab frame, d and m, we can use the Hamilton 
operator for relativistic particles with electric and magnetic dipole moments. For spin-1/2 
particles, it has been derived in Ref. [16]. The Hamilton operator for spin-1 particles is 
similar, because it should be consistent with the corresponding equation of spin motion 
(modified T-BMT equation) which is valid for any spin. 

If we neglect the normal magnetic moment fiQ = eS/m which is small for nuclei, the 
relations between the electromagnetic moments in two frames are given by 

d = d' - ^jP(d' -l3)-(3xm', m = m' - -J—p(m' -f3)+f3xd'. (17) 

When d = el, the relation between d and d! arises from the Lorentz transformation of the 
length of electric dipole, Z, because the charge e is a relativistic invariant. Relations (17) 
remain valid for induced electromagnetic moments. 

As a result, the correction to the Hamilton operator in the rotating frame takes the form 

V = -^- (d' -E' + m'- B') = -— (S ■ E'f - ^{S ■ B'f. (18) 

27 7 7 
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Eq. (13) is an approximate version of Eq. (18). 

The equation of oscillatory motion of the particle has the form 

< 19 > 

The quantity E in Eq. (19) is the electric field met by the particle. As a result of the 
coherent beam oscillations, this field oscillates in the particle's rest frame. The angular 
frequency of velocity modulation uo significantly differs from that of the resonator (see Refs. 
[6, 8]) but should be very close to the angular frequency of spin rotation lo . The latter 
quantity is almost equal to the vertical component of u) a , because other components of this 
vector are relatively small: 

oo, = (u a ) z = B z . (20) 

For the deuteron, lo > 0. 

The modulation of normalized velocity can be given by (see Refs. [6, 8]) 

(3 = f = p + A/3 • cos (cot + cp) e4> , (21) 



where 



a Po \Jm 2 + pl 

A) = ; 70 = 



\Jm? + pi m 
Owing to this modulation, the radial electric field in the particle's rest frame has the oscil- 
latory part. The effect of the modulation on the BVP is described by the last term in Eq. 
(6) proportional to f3xB. 

We can perform the calculation of electric field acting on the particle to within first-order 
terms in A(3 . The particle momentum is defined by the equation 

P= Tyf^ff =Po + 7o mA A) - cos (cot + ip)e 4> . (22) 

According to the result of differentiation of Eq. (22) on time, Eq. (19) takes the form 

E = -E sin (cot + <p )e 4> , (23) 

where 

Eo = 2l!^ AA . ( 24 ) 

e 

Eq. (18) can be transformed to the form 

V = -^(^B-^Sp + E^f - (3tiB 2 z S 2 z . (25) 



An estimate of two terms in the formula for the effective electric field E' [see Eq. (14)] 
shows that the term proportional to the magnetic field B z is much bigger for the deuteron. 
To simplify the calculation, we neglect the effect of the longitudinal electric field and use 
the approximation 

V = - 1 B 2 z (a T [3 2 S 2 p + [3 T S 2 ). (26) 
The quantities 7 and /3 2 7 can be expanded in series of A/3 : 

7 = 7o + A)7o • AA, cos (ut + y?) + 1(1 + 3/3 2 7o)7o(AA)) 2 {1 + cos [2(ut + <p)]} , 

I3 2 l = Pfao + (2 + (3fa 2 )(3 lo ■ Aft cos {oot + y?) (27) 
+ \{2 + 5/3 2 7o 2 + 3/3 4 7o 4 )7o(A/5o) 2 {1 + cos [2{uot + <p)}} . 

Eqs. (26), (27) define the corrections to the Hamilton operator for the tensor polarizabil- 
ities of deuteron. 



V. SOLUTION OF MATRIX HAMILTON EQUATION 



Nonzero matrix elements of the spin operators contained by Eq. (26) are 



(Si 



111 



(Sp)i -1 - (^)-i,i - ('S'p)-i-i - o' ( S p)oo - 1) (^)ii - (-S'z)-i-i - 1- 



where 



Therefore, matrix Hamiltonian (5) takes the form 

( E + u + A + B A 
h= /•;„ • 1A 

A E -u + A + B) 



\ 



A = a + ai cos {ujt + (/?) + a 2 cos [2(ut + <£>)], 
£> = b + cos (uot + (p) + b 2 cos [2(atf + </?)], 



a 



--a T ^7o 



/5 2 + i (2 + 5^ 2 7 2 + 3/5 4 7o 4 )(A/5o) 2 



ax = --a T 5 2 (2 + /3 2 7o 2 )A,7o • AA>, 
« 2 = -i« TJ B 2 (2 + 5/3 2 7o 2 + 3/3 4 7o 4 ) 7o (A/3 ) 2 , 
fl + ^(l + 3/3 2 7 2 ) 7o 2 (A/3 ) 2 



b 2 = --3 T B 2 Z {1 + 3/5 2 7o 2 ) 7o 3 (A^ ) 2 , 



(28) 



(29) 



(30) 
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and Eq is the zero energy level. 

In Hamiltonian (29), the EDM effect is not taken into account. 

We consider the spin dynamics near a resonance. At the first stage, it is useful to pass 
on to new amplitudes (see Ref. [11]). This transformation brings real parts of diagonal 
elements of the matrix Hamiltonian to zero. However, it does not nullify the imagine parts 
of diagonal elements for unstable particles (see Ref. [17]). Evidently, Hamiltonian (29) is 
real and the new amplitudes are equal to 



7i(t) = exp Y 

7o(0 = exp\i 
7_i (t) =exp\i 



j . . fli + h a 2 + b 2 , 
kit + fit) + — g(t) 



uj 2uj 

M + — f(t) + ^g(t) 

, , , ai + h * a 2 + b 2 ( 
k-it + fit) + — g(t) 



Ci(t), 
C (t), 

C-i(t), 



(31) 



uj 2uj 

ki = E + u + a + b , k = E + 2a , = E - u + a + b , 

f(t) = sin (tut + if) — sin (if), g(t) = sin [2(ut + if)] — sin (2if). 
Dynamics of these amplitudes does not depend on the tensor magnetic polarizability and 
is given by 



d'ji 

i—T- = Aexp {2iuJot)7-i 

(JjV 

# = o 

dt 

i — — = -4 exp (— 2tuj t) , fi 



(32) 



dt 



M + — f{t) + -g{t) 

UJ UJ 



Co(0). 



(33) 



Eqs. (31), (32) result in 

C (t) = exp 

Zero component of spin is not mixed with other components. 

At the second stage, we can average over much longer time than the oscillation period 
[11]. The relation 

cos ((t + rj) = i {exp [i((t + rj)} + exp [— i(C* + v)}} 

can be used. There are two resonant frequencies, uj = uj and uj = 2uj q . First of them 
corresponds to the resonance condition in the deuteron EDM experiment [6, 7]. We will 
consider this case first. 

When uj uj q , averaging over time results in 



.dji 

1 dt 
.OT-i 

dt 



a 2 

~2 

a 2 



7_i exp {2i[(uj - uj)t -if]} 
7i exp {-2i[(uj - uj)t -if]} 



(34) 
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At the third stage, we can use the following transformation: 

£>i(f) = exp [-i(u - w)*]7i(*), 
= exp [i(a;o - u;)i]7i(i). 
Transformed equation (34) can be written in the matrix form: 



dt 



( 



UJq — UJ 



a 2 



exp (— 2iy?) 



a 2 



D = 



£>x(f) 
£>_x(t) 



Eq. (36) can be analytically solved. The solution of it has the form 
D i(t) = \ U} ' + ^°~ U} D 1 (0) + ^-p_M exp {-iJt) 
+ — VAO) - —V-i(O) exp(iu't), 



D-i(t) 
+ 



2uj' y ' 2uj' 

|> l(0 ) + ^^) fl . l(0) 



£ > l(0 ) + ^±^ O _ l( 0) 



2w 



exp (—iuj't) 
exp (iu/t) 



or 



D x [t) = cos(cj't) -i^^sin(cj't) IVO) - i— sin (w'*)I>-i(0), 

a/ J uj' 

r 

£>_i(t) = -i— sin (w'*)Di(0) + cos (u't) + i^Z^ s i n (^) D^(0), 
uj L uj 



where 



UJ 



{us - uj) 2 + ££*, S = y ex P (-2*v ? )- 



The angular frequency of spin oscillation is equal to 2uj' 
The initial spin amplitudes take the form 

CLi + b 



Ci(t) = exp 
C-i(t) = exp 



(E + uj + a + bo)t + 
(E - uj + a + &o)* + 



i . a 2 + b 2 , , 



>£>-i(*), 



2cj 

d(0) = ^(0), C_ 1 (0) = J D_ 1 (0). 
The resonance at the doubled frequency ss 2uj can be investigated in a similar 
The evolution of the spin amplitudes is given by 



Ci(t) = exp < —% 



(£ + - + a + bo)t + -f(t) + -^r-^g(t) 

Z UJ ZUJ 

C (t) = exp {-i \(E + 2a )t + ^f(t) + -g(t)] }c (0), 
I L uj uj J J 



C_i(i) = exp < — % 



(£ - 77 + a + bo)t + - -f(t) + -^-^g(t) 

Z UJ ZUJ 

C 1 (0) = D 1 (0), C_ 1 (0) = J D_ 1 (0), 



>£>-i(*), 
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where 



£>x(f) = (cos — - sm _ j Dl(0 ) _ sm — £Li(0), 

£>-i(*) = -i— sin ^T^i(O) + cos — + sin — LL^O), ( 42 ) 

£' = y ex p(-^)> 

and the angular frequency of spin oscillation is equal to 



uj" = J (2uj -uj) 2 + &'£'*. (43) 



VI. SPIN DYNAMICS CAUSED BY TENSOR POLARIZABILITIES OF 
DEUTERON 

For spin-1 particles, three components of polarization vector and related components of 
polarization tensor are defined by 

p — "^7^v° lLy + U l°0 + ^0^-1 + ^0°-lJ' 

3 1 

^4>4> — 1 1 1 ~~ °o^oJ ~~ 2' 

P^ = <?(C7 1 C* 1 -C?C7_ 1 ). 

The horizontal components, P p and P^, do not give necessary information about the 
investigated effect because they undergo fast oscillations caused by the g— 2 spin precession. 
The change of the vertical component, P z , is a relatively slow process. 

The quantity P z does not depend on C . Since C\C{ = DiD{, C-\C*_ X = D-iD*_ ± , the 
BVP is caused by the tensor electric polarizability and is not affected by the tensor magnetic 
one. However, this conclusion is not valid if the deuteron possesses the EDM. In this case, 
the tensor magnetic polarizability leads to splitting of resonance frequency [1, 2, 3]. 

When uj m uj , the evolution of the vertical component of polarization vector is expressed 
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by 



1 - -|(1 -cos(2u/f)) 



\ ^[Pp P (0) - P H (0)\ cos (2<p) (1 - cos (2u't)) - sin (2<p) sin (2u't) 

+P p4 ,(0) U ° ~ - sin (2<p) (1 - cos (2u't)) + cos (2<p) sin (2u't) } , £ = ^, 



(45) 



where the quantities a 2 and a;' are defined by Eqs. (30) and (39), respectively. 

For a vector-polarized deuteron beam, the related components of polarization vector and 
polarization tensor have the form 



P z = cos(9), P, 



pp 



3sin 2 (fl) cos 2 (ip) 



P _! 



3 sin 2 (9) sin 2 (ip) - l] , P p<t> = - sin 2 (0) sin (2ip) , 



(46) 



where 9 and ip are the spherical angles defining the spin direction in the rotating frame. The 
azimuth ip = characterizes the spin directed radially outward. 

The projection of deuteron spin onto the preferential direction can be equal to zero. The 
beam possessing such a polarization is tensor polarized. The vector polarization of this beam 
is zero. The components of polarization vector and polarization tensor take the form 



P p = P 4> = P z = 0, P pp = -3 sin 2 (9) cos 2 (V>) + 1, 
P H = -3 sin 2 (9) sin 2 {$) + 1, P p4> = ~ sin 2 (0) sin (2^), 



(47) 



where 9 and ip are the spherical angles stated above. When the polarization of deuteron 
beam is vector, Eqs. (45), (46) result in 

E 2 



PM 



sin 2 (9) ( — - ^ cos [2(^ - <p)] [1 - cos (2c/t)] + sin [2(^ - </?)] sin (2«/t)} . 
2uj I w J 



1 - 



/2 



1 -cos(2u/t)) 



cos (0) 



(48) 



If the resonance is perfect, the initial beam polarization is horizontal [P z (0) = 0], and 
uj't « 0, Eq. (45) takes the form 



P z (t) = -a 2 t 



P P M cos (2^) - Ppp(0) - P ^ (0) sin (2^) 



(49) 



In this case, the vertical spin component grows linearly with time and 



dP z 2 



— ~a>2 
dt 3 



P p ,(0) cos (2^) - PpM - P(t>M sin (2 V ) 



(50) 
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For the resonance at the doubled frequency uo 2u; , the evolution of the vertical com- 
ponent of polarization vector is given by 



P*{t) = 



If 

+^ \ [Ppp(O) - P*(0)] 



4£ 



/2 



3cu' 
+2P*(0) 



2cuq — c<j 



2a;o — cu 
77' 



u2 



?(l-cos(u; ff f)) 



cos (</?) (1 — cos ((*/'£)) — sin (</?) sin (w"t) 



(51) 



a;" 



sin ((f) (1 — cos (u/'i)) + cos (</?) sin (u/'t) 



ai 



where the quantities a\ and u;" are defined by Eqs. (30) and (43), respectively. 

If the resonance is perfect, the initial beam polarization is horizontal, and uj"t « 0, Eq. 
(51) takes the form 



P z (t) = - ai t 



Pp^O) cos - — ^ — sin (</?) 



(52) 



Eqs. (30), (49), and (52) show that resonance frequency doubling leads to a dramatic 
amplification of the Baryshevsky effect. When the frequency is doubled, the EDM effect 
becomes nonresonant. In this case, it does not influence the spin dynamics. 



VII. MEASUREMENT OF TENSOR ELECTRIC POL ARIZ ABILITY OF 
DEUTERON IN STORAGE RING EXPERIMENTS 

To discover the Baryshevsky effect, it is necessary to stimulate the BVP conditioned by 
the tensor electric polarizability of deuteron and to avoid a similar effect caused by the 
magnetic moment. It is known that the magnetic resonance takes place when the particle 
placed in a uniform vertical magnetic field is also affected by a horizontal magnetic field 
oscillating at a frequency close to the frequency of spin rotation (see, e.g., Ref. [20]). The 
magnetic resonance results in the BVP for a horizontally polarized beam. 

Evidently, the magnetic resonance cannot take place when the electric field is longitu- 
dinal, because nothing but the oscillating electric field appears in the particle's rest frame. 
Since the frequencies of betatron oscillations are chosen to be far from resonances, these 
oscillations cannot lead to the resonance effect. However, the resonance is caused by the 
tensor electric polarizability of deuteron. The electric field in the particle's rest frame pos- 
sesses the oscillating longitudinal component and the radial one E' caused by the Lorentz 
transformation of the vertical magnetic field. The latter component has a resonance part 
because of the modulation of the particle velocity. 
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To measure the effect, some resonators (rf cavities) should be used. The electric field 
in a resonator is generated along the central line and the magnetic field is orthogonally 
directed [21]. The magnetic field along the central line is equal to zero. If the rf cavities 
are perfectly placed and longitudinally directed, the magnetic field cannot stimulate any 
resonance effect. In this case, the observed BVP corresponds to the definite value of the 
tensor electric polarizability. 

However, both a displacement and an angular deviation of the center line of the rf cavities 
away from an average particle trajectory lead to a similar behavior of spin imitating the 
Baryshevsky effect. As a result, they create systematical errors in measurement of the tensor 
electric polarizability. Most of these errors are not in resonance with the spin precession 
in the horizontal plane. Therefore, they create background and result in fast oscillations of 
the vertical component of polarization vector (see Refs. [6, 7, 8, 22]). Besides this effect, 
the systematical error can be caused by a radial magnetic field in the particle's rest frame 
oscillating at the resonant frequency. In the deuteron EDM experiment, a similar error 
will be eliminated by alternately producing two sub-beams with different betatron tunes 
[6, 8, 22]. In the deuteron tensor-electric-polarizability one, the resonant radial magnetic 
field in the particle's rest frame is much less important when a tensor polarized beam is 
used (see below). We calculate only effects of resonant fields on the BVP in ideal conditions 
and disregard systematical errors. Thus, we take into consideration the constant vertical 
magnetic field and the oscillating longitudinal electric one. 

The measurement of the tensor electric polarizability of deuteron in a storage ring needs 
the field configuration similar to that proposed for the deuteron EDM experiment [6, 7, 8]. 
However, the resonance frequency should be doubled (u ~ 2cj )- Resonance frequency 
doubling cannot be implemented in the designed EDM ring. In this ring, the eigenfrequency 
of free synchrotron oscillations must be chosen close to the g— 2 frequency, u a , and the 
resonance effect is created by the beatings between two rf frequencies [8]. Therefore, the 
measurement of the tensor electric polarizability of deuteron needs another ring or at least 
rf cavities different from that developed for the deuteron EDM experiment. 

However, the Baryshevsky effect caused by the tensor electric polarizability of deuteron 
should be taken into account when performing the deuteron EDM experiment [1, 2, 3]. This 
effect results in the similar BVP and can imitate the presence of the deuteron EDM of order 
of d ~ 1CT 29 e-cm. An attainment of such an accuracy is the goal of the storage ring EDM 
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experiment [6, 7, 8]. 

The EDM-dependent evolution of deuteron spin in this experiment has been calculated 
in detail in Ref. [23]. The dynamics of the vertical component of polarization vector is given 
by 

p(EDM ){t) = g j^L^ cos -<p)[l- cos (n't)} + sin (V - <p) sin (Q>t)} , (53) 

where 

n' = \Q'\ = - ujf + Sf, (54) 

£S = -\dB z ■ A/3 f 1 + ^t) , (55) 

and the azimuth ip defines the direction of spin at zero time. The initial polarization is 
supposed to be horizontal. 
When Q't « 1, 

p(EDM) =£ . tsin (^_ ip ) 

= ~dB z A(3 (\ + ^L] tsm(^ - <p). 

Z \ LO J 

We can evaluate the expected sensitivity in the measurement of the tensor electric polar- 
izability of deuteron with the comparison of Eqs. (51) and (53) and the use of Eq. (46) and 
the sensitivity of the deuteron EDM experiment estimated in Ref. [8]. For the deuteron, 
a = a d = -0.14299. The sensitivity to the EDM of d — 1 x 1(T 29 e-cm corresponds to the 
accuracy of Sax = 1.2 x 10~ 43 cm 3 when cu w 2cu and Eqs. (51),(53)-(55) are used. This 
estimate is based on the values of 70 = 1.28, (3q = 0.625, At> = 3.5 x 10 6 m/s, and B z = 3 
T [8]. There are three independent theoretical predictions for the value of the tensor electric 
polarizability of deuteron, namely ax = —6.2 x 10~ cm [24], -6.8 x 10~ 41 cm 3 [25], and 
3.2 x 10~ 41 cm 3 [26]. Two first values are very close to each other but they do not agree 
with the last result. 

In all probability, the best sensitivity in the measurement of ax can be achieved with 
the use of a tensor polarized deuteron beam. The initial preferential direction of deuteron 
polarization should be horizontal. When the vector polarization of such a beam is zero, any 
spin rotation does not occur. In this case, there are no related systematical errors caused 
by the radial magnetic field and some other reasons. In the general case, such systematical 
errors are proportional to a residual vector polarization of the beam. This advantage leads 
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to a sufficient increase in experimental accuracy. When u ~ 2ujq, the equation describing 
the evolution of beam polarization takes the form 

P z (t) = sin 2 cos (2V - *>) [1 - cos (aA)] 

[ a; (57) 

+ sin (2^ — 99) sin (u"t) >. 

In this case, the preliminary estimate of experimental accuracy is Sax ~ 1CT 45 -r- 10~ 44 
cm 3 . 

When = 7r/2, the natural choice of phase 

^ , , 7T 

V9 = 2^ ± - 

brings Eq. (57) to the form 

P z (t) = ±^sm(uj"t). (58) 

Other possibilities, {p = 2ip and = 2^ ± 7r, lead to the equation 

4£'(2cj - u) 9 /V'f 

^ — sm u 



P.(0 = ^ V rSi ^ T (59) 



The dependence of P z (t) on time becomes quadratic when uu"t <C 1. Therefore, these 
possibilities are less useful. However, they can be used for checking the result. 

The deuteron tensor-electric-polarizability experiment essentially differs from the 
deuteron EDM one by a nonnecessity of careful checking systematical errors caused by 
the horizontal magnetic field in the particle's rest frame. The use of a tensor polarized beam 
makes it possible to avoid any spin rotations and to cancel all related systematical errors. 
The BVP of such a beam is defined only by the tensor electric polarizability of deuteron. 
It is a great advantage because the elimination of similar systematical errors is one of main 
problems for the deuteron EDM experiment [8, 22]. A residual vector polarization of the 
beam together with a resonant magnetic field in the particle's rest frame can result in a false 
signal. However, a necessary correction into the BVP can be made with a longitudinally 
vector-polarized beam. It is important that the BVP caused by any systematical error is 
overturned and the BVP defined by the tensor electric polarizability remains the same when 
reversing the polarization of this beam [see Eq. (48)]. This property brings a easy differ- 
entiation between the Baryshevsky effect and false signals for vector-polarized beams. In 
all probability, the deuteron tensor-electric-polarizability experiment can be made at one of 
existing rings. 
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VIII. DIFFERENTIATION OF EFFECTS OF EDM AND TENSOR ELECTRIC 
POL ARIZ ABILITY IN THE DEUTERON EDM EXPERIMENT 



Eqs. (45) and (53)-(55) describing the effects of the tensor electric polarizability and 
the EDM on the spin dynamics in the EDM experiment essentially differ. Therefore, these 
effects can be differentiated. 

When the initial polarization of deuteron beam is horizontal, Eq. (48) takes the form 

pitensor)^ = f ^ ^ _ ^ ^ _ ^ + ^ ^ _ ^ ^ _ ^ 

2cu L lo J 
For the EDM experiment, the choice of phase 

<p = i> ± | 

is necessary. This choice results in 

p(EDM) {t) = ± g gin ^ p{tensor) {t) = ± ^ - Uj) ^ ^ ^ 

Since the quantities Sq,S are very small, Q' ~ a;' ~ |cj — <^|- When we analyze only 
the ratio of amplitudes in Eq. (61) which is approximately equal to 2£q/£ , we obtain 
that the values of «t found in Refs. [24, 25, 26] correspond to the false EDM moments of 
\d\ — 3 x 1CT 29 , 3 x 1CT 29 , and 2 x 1CT 29 e-cm, respectively. However, the EDM contribution 
to P z grows linearly with time when Q't <C 1, while the tensor-electric-polarizability contri- 
bution is negligible in this case. Therefore, keeping the frequency and phase of the coherent 
longitudinal oscillations almost equal to the frequency and phase of the spin rotations makes 
it possible to cancel the effect of the tensor electric polarizability in the framework of the 
deuteron EDM experiment. The same conclusion was recently drawn in Ref. [27]. Never- 
theless, the tensor electric polarizability of deuteron should be taken into account. To check 
the possible existence of the deuteron EDM, one can also use other possibilities of separating 
the EDM and Baryshevsky effects listed below. 

1) The spin dynamics caused by first-order interactions (including the EDM effect) and 
second-order interactions (including the Baryshevsky effect) is defined by the operator equa- 
tions of spin motion 

y = AfixS (62) 

and 

—j^ = PijkSjSk, (63) 
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respectively. Therefore, the EDM effect reverses the sign when the beam polarization is 
reversed while the sign of the Baryshevsky effect remains unchanged. 

2) Since both the EDM and Baryshevsky effects depend on the difference ift — tp, reversing 
the beam polarization (/0 — > ip + 7r) is equivalent to the transition to the opposite phase 
(</?—></? + 7r) . Naturally, such a transition is technically simpler. 

If two measurements fulfilled according to points 1) or 2) give the values P z \ and P z2 for 
the BVP, the EDM and Baryshevsky effects are characterized by the values (P z i — P z2 )/2 
and (P z i + P z2 )/2, respectively. 

3) In the particle rest frame, the EDM and Baryshevsky effects are linear and quadratic 
in the electric field, respectively. The experimental dependence can be determined with 
changing the amplitude of the field in the resonators. 

4) The frequency of BVP caused by the Baryshevsky effect is approximately twice as 
large as that conditioned by the EDM. 

5) The use of tensor polarized deuteron beam even at the angular frequency uj ~ uj 
cancels the EDM effect and main systematical errors. The evolution of beam polarization 
is given by 

pitensor)^ = (^J^ CQS fity _ ^)] [1 _ COS (2u't)\ + SU1 [2(^ - if)] sin (2w't)} , 

2oj I w J 

(64) 

if the initial polarization is defined by Eq. (47). 

Thus, the EDM and Baryshevsky effects can be effectively differentiated and the latter 
effect can be cancelled in the framework of the deuteron EDM experiment. 

IX. DISCUSSION AND SUMMARY 

The spin dynamics conditioned by the tensor electric polarizability of deuteron has been 
calculated for the first time in Refs. [f, 2, 3]. To compare our results with those obtained 
in [f , 2, 3], it is helpful to introduce the effective field defined by 

E 2 eff = P 2 JBI (65) 

In Refs. [f, 2, 3], the effect has been described to within first-order terms in (3. In this 
approximation, the squared effective field is equal to 

E 2 eff = (Ei%f + 2B% ■ AA) cos (ut + <p) + ^B 2 Z (A(3 ) 2 cos [2{ut + <p)\. (66) 

20 



The evolution of polarization vector is given by Eqs. (24), (29) in Ref. [2] and Eqs. 
(44), (49) in Ref. [3]. The final equation has the form 



dP 1 

2 = — Att T cos (2Q f t + 2(p f ) 



dt 2 
where 



P^(0) cos (2Sti) - Ppp(0) P ^ (0) sin (2Qt) 



(67) 



A q t = -^a T B 2 z (A(3 ) 2 , (68) 



n corresponds to our designation ujq, and (pf = (p + ix/2. The spin rotation is supposed to be 
clockwise in Refs. [1, 2, 3] and counter-clockwise (uj > 0) in the present work. Therefore, 
Qf ps — Q and averaging Eq. (67) over time with allowance for Eq. (68) results in 

dP 1 

_£ = __ aTjB 2 (AA)) 2 {2 p^( ) cos (2^) - [P P M - P H (0)\ sm (2^)} . (69) 

This equation fully agrees with Eq. (50). 

Agreement of the results obtained in Refs. [2, 3] and in the present work confirms their 
validity. The method used in Refs. [2, 3] is less convenient to calculate the spin dynamics 
in oscillatory external fields than in static ones. In the theory of magnetic resonance, the 
transition to a rotating frame is commonly used [20]. Evidently, the transition to the rotating 
frame is necessary in order to determine the spin dynamics with the method developed in 
Refs. [2, 3] when u ^ uj (Q ^ \flf\). 

The calculated effect of the BVP caused by the tensor electric polarizability of deuteron 
is an exciting example of new spin physics brought by tensor interactions. In the considered 
case, the deuteron spin is governed by the electromagnetic interaction. The similar effect 
investigated in works by V. Baryshevsky [28, 29] is affected by the strong interaction of 
deuteron with nuclear matter. These effects stimulated by tensor interactions result in the 
transformation of tensor polarization into the vector one and the other way round. 

Eqs. (37) shows that the spin-up state is converted into the spin-down state and the 
other way round. This property is caused by the nondiagonal terms in Hamiltonian (36). 
As a result, the vertical component of polarization vector oscillates. This phenomenon is 
similar to light birefringence in crystals [28, 29]. 

A similar behavior of spin takes place at the NMR when a nucleus is placed into a 
resonant horizontal magnetic field. As is well known, the NMR also consists in an oscillation 
of P z . However, there is an essential difference between two effects. The Baryshevsky effect 
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exists even when the beam is tensor polarized, while the NMR does not change the beam 
polarization in this case. 

The calculation shows that the Baryshevsky effect can be observed in storage rings. 
Performing the measurements with the use of resonance u ~ 2ujq offers an opportunity to 
measure the deuteron's tensor electric polarizability with the accuracy of 10 -45 -j- 1CT 44 cm 3 
(1(T 6 ^ 1(T 5 fm 3 ). 

It is also possible to use low-energy deuterons in a Penning trap. 

The problem of influence of the tensor electric polarizability on spin dynamics in the 
deuteron EDM experiment in storage rings has been investigated. The EDM and Bary- 
shevsky effects can be effectively differentiated and the latter effect can be cancelled in the 
framework of this experiment. 

In the present work, general formulae describing the BVP conditioned by the tensor 
electric polarizability have been derived. Calculated formulae agree with the previous results 
[1, 2, 3] obtained in a more particular case. The method based on the use of Hamiltonian 
approach and spin wave functions happens to be very convenient for the investigation of the 
effect. 
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